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Abstract 

We study the upper hmits on the mass of the hghtest and sec- 
ond Ughtest CP even Higgs bosons in the (M+1)SSM, the MSSM 
extended by a gauge singlet. The dominant two loop contributions 
to the effective potential are included, which reduce the Higgs masses 
by ~ 10 GeV. Since the coupling R of the lightest Higgs scalar to 
gauge bosons can be small, we study in detail the relations between 
the masses and couplings of both lightest scalars. We present upper 
bounds on the mass of a 'strongly' coupled Higgs (i?>l/2)asa func- 
tion of lower experimental limits on the mass of a 'weakly' coupled 
Higgs {R < 1/2). With the help of these results, the whole parameter 
space of the model can be covered by Higgs boson searches. 
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1 Introduction 



Curiously enough, the most model independent prediction of supersymmetric 
extensions of the Standard Model concerns a 'standard' particle: the mass rrih 
of the (lightest CP even) Higgs boson. Within the minimal supersymmetric 
extension of the Standard Model (MSSM) its mass is bounded, at tree level, 
by 

ml < Ml cos'^2P (1.1) 

where tan/3 = hi/h2 {Hi couples to up-type quarks in our convention). It has 
been realized already some time ago that loop corrections weaken this upper 
bound m . These loop corrections depend on the top quark Yukawa coupling 
ht and the soft susy breaking parameters as the stop masses of 0{Msusy)- At 
the one loop level, given the present experimental errors on the top mass rrit 
and assuming Mgusy ^ 1 TeV, the upper limit on is < 140 GeV. Also 
two loop corrections to rrih have been considered in the MSSM these 
have the tendency to lower the upper bound on by ~ 10 GeV. 

The subject of the present paper is the next-to-minimal supersymmetric 
extension of the Standard Model ((M+1)SSM) |^-|TB[ where a gauge singlet 



superfield S is added to the Higgs sector. It allows to omit the so-called 
/i term fiHiH2 in the superpotential of the MSSM, and to replace it by a 
Yukawa coupling (plus a singlet self coupling): 



W = XSHiH2 + -S^ + ... (1.2) 
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The superpotential (|1.2|) is thus scale invariant, and the electroweak scale 
appears only through the susy breaking terms. 

In view of ongoing Higgs searches at LEP2 [|l^-|l6| and, in the near future. 



at Tevatron Run II ^7j, it is important to check the model dependence of 
bounds on the Higgs mass. In the (M+1)SSM, the upper bound on the mass 
nil of the lightest CP even HiggsQ differs from the one of the MSSM already 
at tree level: now we have 

/ 2A^ \ 

< M| cos^2/5 + — 2 sin22/5 (1.3) 

V 9i + 92 J 

where gi and g2 denote the U{1)y and the SU(2)l gauge couplings. Note 
that, for A < .53, mi is still bounded by Mz at tree level. Large values of A, 



^As there are three CP even Higgs states in the (M+1)SSM, we denote them as Si with 
masses m^, i=1..3, in increasing order. 
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A > .7, are in any case prohibited, if one requires the absence of a Landau 
singularity for A below the GUT scale P,^]. 

Loop corrections to mi have also been considered in the (M+1)SSM 
0. Given rrit and assuming again Mgusy < 1 TeV, the upper limit on mi 
at one loop is then ~ 150 GeV. Within the constrained (M+1)SSM (the 
C(M+1)SSM), where universal soft susy breaking terms at the GUT scale 
are assumed [[7|-[lO[|, A is always below ~ .3, and the upper limit on mi reduces 
to the one of the MSSM (at one loop) of ~ 140 GeV. Two loop corrections 



in the (M+1)SSM have recently been considered in ITS 



Within the (M+1)SSM this is, however, not the end of the story: It is 
well known 0,[iy,[ll]] that now the lightest Higgs scalar 5*1 can be dominantly 
a gauge singlet state. In this case it decouples from the gauge bosons and 
becomes invisible in Higgs production processes, and the lightest visible Higgs 
boson is then actually the second lightest one 5*2. Fortunately, under these 



circumstances 5*2 cannot be too heavy |]^, |T0|, |TT| : In the extreme case of a pure 
singlet lightest Higgs, the mass m2 of the next-to-lightest Higgs scalar is again 
below the upper limit designed originally for mi. In general, however, mixed 
scenarios can be realized, with a weakly coupled (but not pure singlet) lightest 
Higgs, and a second lightest Higgs above the previous mi limits. Although 
analyses of the Higgs sector including these scenarios in the (M+1)SSM have 
been presented before [|rT| we find that these should be improved: First, 
experimental errors on the top quark pole mass mf°'^ = 173.8 ±5.2 GeV |18] 



have been reduced considerably, leading to stronger constraints on the top 
quark Yukawa coupling ht which determines to a large extent the radiative 
corrections to mi^2- Second, at least the dominant two loop corrections to the 
effective potential should be taken into account, since they are not necessarily 
negligible. The purpose of the present paper is thus an analysis of the allowed 
masses and couplings to the gauge bosons of the lightest CP even Higgs 
scalars in the (M+1)SSM, including present constraints on and a two 
loop improvement of the Higgs potential. 

In the next section we present our method of obtaining the dominant two 
loop terms in the effective potential, and in section ^ we give the resulting 
upper bound on the lightest Higgs mass. Albeit this upper limit can be ob- 
tained analytically, the mass of the second lightest Higgs in relation to the 
coupling to the gauge bosons requires a numerical analysis. Our methods of 
scanning the parameter space of the model in two different scenarios (con- 
strained and general (M+1)SSM) are presented in section ^. Results on the 
Higgs masses and couplings, and conclusions are presented in section ^ 
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2 Two loop corrections 



In order to obtain the correct upper limit on the Higgs boson mass in the 
presence of soft susy breaking terms, radiative corrections to several terms 
in the effective action have to be considered. Let us first introduce a scale 
Q ~ Msusy, where Mgusy is of the order of the susy breaking terms. Let 
us assume that quantum corrections involving momenta > have been 
evaluated; the resulting effective action re//((5) is then still of the standard 
supersymmetric form plus soft susy breaking terms. Assuming correctly nor- 
malized kinetic terms (after appropriate rescaling of the fields), the Q de- 
pendence of the parameters in FeffiQ) is given by the supersymmetric /3 
functions (valid up to a possible GUT scale Mqut)- 

Often one is interested in relating the parameters in TeffiQ) to more 
fundamental parameters at Mqut- To this end one integrates the supersym- 
metric renormalization group equations between Mqut and Q ~ M^usy to 
one or, if one whishes, to two loop accuracy. Note, however, that the limits 
on the Higgs boson mass depend exclusively on the parameters in TeffiQ) at 
the scale Q ~ Mgusy] the two loop contributions to the effective potential con- 
sidered below serve to specify this dependence more precisely. The accuracy 
to which one has (possibly) related the parameters at the scale Q ~ Mg^gy to 
parameters at a scale Mqut is completely irrelevant for the relation between 
the Higgs boson mass and the parameters at the scale Q ~ Msusy 

One is left with the computation of quantum corrections to Fg// involving 
momenta ^ Q^- Subsequently the quantum corrections to the following 
terms in Fg// will play a role: 

a) Corrections to the kinetic terms of the Higgs bosons. Due to gauge in- 
variancc the same quantum corrections contribute to the kinetic energy 
and to the Higgs-Z boson couplings, which affect the relation between 
the Higgs vevs and M^; 

b) Corrections to the Higgs-top quark Yukawa coupling; 

c) Corrections to the Higgs effective potential. These corrections could, in 
principle, be decomposed into contributions to the Yukawa couplings 
A and /t of eq. (1.2) and the soft terms (these contributions are the 
ones proportional to InQ^ or, at two loop order, In^Q^), and "non- 
supersymmetric" contributions which are Q"^ independent. These latter 
contributions to the effective potential are of the orders (vev)^ with 
n > A and become small in the case of large soft terms compared to 
the vevs. Our results in section 5 are based on the effective potential 
including these contributions (which are not necessarily numerically 
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irrelevant), and there is no need to perform the decomposition of the 
radiative corrections to the effective potential explicitely. 

Let us start with the last item: The Higgs effective potential Vejf can be 
developped in power of h or loops as 

V,jf = V^W + r(^) + V^(2) + . . . . (2.1) 

Within the (M+1)SSM, we are interested in the dependence of Vejf in three 
CP even scalar vevs hi, h2 and s (assuming no CP violation in the Higgs 
sector). The tree level potential V^^^-* is determined by the superpotential 
( |1.2| ) and the standard soft susy breaking terms []5|- p]] . For completeness, 
and in order to fix our conventions, we give here the expression for V^^^: 



V^(o) 



'mjj^hi+ mjj^h2 + mgS - 2\Axhih2S + - kA^^s 



2„4 



8 



2\2 



{hi - hi) 



The one loop corrections to the effective potential are given by 

1 



^(1) 



6471^ 



STrM^ 



In 



(2.2) 



(2.3) 



where we only take top and stop loops into account. The relevant field 
dependent masses are the top quark mass 



mt = hthi 

and the stop mass matrix (in the {T^, Tl) basis) 



ml + ml 



(2.4) 



(2.5) 



where m^, are the stop soft masses and 

At = At-\scoi^ (2.6) 

is the so-called stop mixing. In eq. ( p.5|) we have neglected the electroweak 
D terms which would only give small contributions to the effective potential 
in the relevant region m-p, mq ^ Mz- The masses of the physical eigenstates 
tl, t2 then read 



4 



m 



(2.7) 



with M; 



susy 



{ml + m%) 



and 5 



itlq + 



(2. 



Note that the top Yukawa couphng ht in eq. ( ^.4| ) and below is defined at 
the scale Q, cf. the discussion at the beginning of this section. 

In the case of large susy breaking terms compared to the vevs /ij, V^'^^ can 
be expanded in (even) powers of hi. The terms quadratic in hi will not affect 
the upper bound on the Higgs mass (and can be absorbed into the unknown 
soft parameters m/^-j, and Ax in ( |2.2| )). In the approximation where the 
stop mass splitting 5 is smallQ, the quartic terms read 



yd) 



3hl 
167r2 



(2.9) 



where t = ln 



and Xt = 2 



1 - 



A? 



12(Mi,^ + m?)^ 



(2.10) 
(2.11) 



In our computations, however, we used the full expression (OD for we 
will use the quartic terms ( |2.9|) in the next section only in order to compare 
our two loop result to those of refs. [§,|l^]. 

Next, we consider the dominant two loop corrections. These will be nu- 
merically important only for large susy breaking terms compared to hi, hence 
we will expand again in powers of hi. Since the terms quadratic in hi can 
again be absorbed into the tree level soft terms, we just consider the quartic 
terms, and here only those which are proportional to large couplings: terms 
~ aghf and ~ h^. Finally, we are only interested in leading logs (terms 
quadratic in t). The corresponding expression for V^'^^ can be obtained from 
the explicit two loop calculation of Veff in or, as we have checked ex- 
plicitely, from the requirement that the complete effective potential has to 
satisfy the renormalization group equations also at scales Q < Mgusy, pro- 
vided the non-supersymmetric (3 function for ht is used. One obtains in both 



cases 



^This approximation is well motivated in the C(M+1)SSM where we take universal soft 
terms at the GUT scale. On the other hand, we have checked numerically that, in the 
general (M+1)SSM, the lightest Higgs mass takes its maximal value for 5 ~ 0. 
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V^i? = 3(^) ht{32.a.-h^t\ (2.12) 

Now, we turn to the quantum corrections to the Higgs boson kinetic 
terms. They lead to a wave function renormahzation factor Zh^ in front of 



the DaHiD^^Hi term with, to order h"^ 



1 1 



ZH. = l + 3j^t (2.13) 

Finally, the quantum corrections to the i^i-top quark Yukawa coupling 
ht have to be considered. After an appropriate rescahng of the Hi and top 
quark fields in order to render their kinetic terms properly normalized, these 
quantum corrections lead to an effective coupling htijnt) with, to orders /i^, 

ht{m,) = K{Q) (l + ^ (327ra, - ^ /i^) . (2.14) 



In eqs. ( p.l3|) and ( p.l4|) the large logarithm t is actually given by In 



Q 



where acts as a UV cutoff, cf. the discussion at the beginning of this 
section. In the relevant region Mg^sy ^ f^t the expression ( |2.1CI| ) for t can be 



used here as well. The (running) top quark mass is then given by 

mt{mt) = ht{mt)Z]l\i (2.15) 
and the relation between the pole and running mass, to order a^, reads 

= + (2.16) 

3 Upper bound on the lightest Higgs mass 

In this section we derive an analytic upper bound on the mass of the lightest 
Higgs scalar. First, we summarize our contributions to the effective potential. 
As it is already known, in the (M+1)SSM the upper bound on the lightest 
Higgs mass mi is saturated when its singlet component vanishes 0,|T^,|TT],|TB[. 
One is then only interested in the /ij-dependent part of the effective potential. 
Assuming hi <^ Mgusy, i-e. up to O (/if), one obtains from eqs. ( p.2|) , (^]9|) 
and (CT) 



Veff{hi,h2) = mlhl + mlhl-mlhih2 + ^^{hl-hlf 
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+Hl6^J /^t (^327ra, - - Z.,^ j (3.1) 
with 

— "2 2 I \ 2 2 I J 

= rrijj^ + A s + rad. corrs. , 
= + •^^■^^ + ''^(^d. corrs. , (3-2) 

ml = 2\s{Ax + Ks) + rad. corrs.. 

The radiative corrections in ( p.2|) stem from the contributions to V^^^ and 
V^"^^ quadratic in hi. In the large tan/5 regime (which saturates the upper 
bound on the hghtest Higgs in the MSSM), one is left with only one non- 
singlet light Higgs hi and (|3.1|) simplifies to 

Veffihi) = mlhl + Xhi (3.3) 

with 



A=«±ii + i^? fl^ + A+am (32.a.^lh^)e. (3,4) 



.2 I ^2 ^^^2 /I X / a2 \ 2 

8 ' 167r"2V2"'' ' 7 ' ^Vl67r2/ 2 

(Note that in the large tan/3 regime At = At and no dependence on the 
(M+1)SSM coupling A is left in A.) Now, we can change the variable hi and 
replace it by a variable h[ in terms of which the kinetic term is properly 
normalized, so that we have 



M| = /if. (3.5) 



From eq. ( p.l3|) one finds 



In terms of h'l the effective potential reads 

Veff{h[) = m'^h'^ + A'/if (3.7) 

with 

,2 \ / q!,2 \2 
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Second, recall that ht in the one loop contribution to eq. ( |3.1| ) is given by the 
Yukawa coupling at the scale Q. Hence, we can replace ht{= ht{Q)) in A' by 
htirrit) using eq. ( |2.14D , which allows to relate it directly to the running top 
quark mass. Eq. ( |2.15|) now reads mt{mt) = ht{mt)h[. 

From ( |3.7|) , one obtains the mass of the lightest non-singlet Higgs in 
the case where the singlet decouples (and in the large tan/3 regime) 



1 oPK// 
2 



dh'{ 



12 



(3.9) 



This is just the correct running Higgs mass, but does not include the pole 
mass corrections, which involve no large logarithms and which we will neglect 
throughout this paper. Using ( |3.5|) and expanding A' to the appropriate 
powers of t, the expression for m\ becomes^ 



< = Af|ll-^tj (3.10) 

(i .7, + t + ^ (5 A? - 32.„.) (.Y. + t)t) 

which agrees with the MSSM result in |Q. (Note, however, that the coefficient 
of the term ~ Xtt on the right hand side of (|3.10|) is not necessarily correct, 
since we would obtain terms of the same order if we would take into account 
simple logarithms in the two loop correction V^'^'^ to the potential.) 

The same procedure can be applied for general values of tan/3. Then, one 
has to consider the 2x2 mass matrix ^(dh^dh^Veff) , i,j = 1,2 where the hi are 
properly normalized. Its smallest eigenvalue gives the following upper bound 
on the mass mi of the lightest Higgs boson for arbitrary mixings among the 
3 states {hi,h2,s) (which can be saturated if the lightest Higgs boson 
has a vanishing singlet component) 



/ 2A^ ^ 

mj < M| f cos22/3 + -^r:^ sin22/5 



(3.11) 



+^^^m?(m,)sin2/3 (Ixt + t 



47r2 



167r2 



^ h^t - 327ra, ) (X* + t)t 



The only difference between the MSSM bound Q and ( p.ll|) is the 'tree 
level' term ~ A^ sin^2/?. This term is important for moderate values of tan/?. 



^In eq. ( 3.1(]| ) and below in eq. ( 3.11 ) we omit the argument of ht wherever its choice 
corresponds to a higher order effect. 
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Hence, the maximum of the hghtest Higgs mass in the (M+1)SSM is not 
obtained for large tan/5 as in the MSSM, but rather for moderate tan/5 (as 
confirmed by our numerical analysis, cf. section ^). On the other hand, the 
radiative corrections are identical in the (M+1)SSM and in the MSSM. In 
particular, the linear dependence in Xf is the same in both models. Hence, 
from eq. (|2.11 ), the upper bound on ml is maximized for Xt = 6 (corre- 
sponding to At = ^/6Msusy, the 'maximal mixing' case), and minimized for 
Xt = (corresponding to At = 0, the 'no mixing' case). 



4 Parametrization of the (M+1)SSM 

Eq. ( p. Ill ) gives an upper bound on the lightest Higgs mass mi regardless of 
its coupling to the gauge bosons. In the extreme case of a pure singlet lightest 
Higgs, the next-to-lightest Higgs is non-singlet and the upper bound (|3.11|) 
actually applies to On the other hand, it can occur that the lightest 
Higgs is weakly coupled to gauge bosons (without being a pure singlet) and 
m2 is above the limit (|3.11| ). This case requires a numerical analysis, which 
will be performed in the next section. First, let us present our methods of 
scanning the parameter space of the (M+1)SSM. 

Not counting the known gauge couplings, the parameters of the model 

are 

A , K, ht, Ax, A^, At, mjj^ , m^^ , m| , , (4.1) 

where ht is eventually fixed by the top mass and an overall scale of the 
dimensionful parameters by the Z mass. Now, let us see how to handle this 
high dimensional parameter space in two different scenarios. 



4.1 Constrained (M+1)SSM 

In the C(M+1)SSM the soft terms are assumed universal at the GUT scale, 
the global minimum of the effective potential has to be the global minimum 
and present experimental constraints on the sparticle and Higgs masses are 
applied. The free parameters can be chosen as the GUT scale dimensionless 
parameters 



Ao m, 



Ao, no, hto, ^ (4.2) 

where Aq, M1/2 and ttlq are the universal trilinar coupling, gaugino mass and 
scalar mass respectively. In order to scan the 5-dimensional parameter space 
of the C(M+1)SSM, we proceed as in refs. 
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First, we scan over the GUT scale parameters ( |4.2| ) and integrate numer- 
ically the renormalization group equations ||^ down to the susy scale in each 
case. 

Then, we minimize the complete two loop effective potential in order to 
obtain the Higgs vevs hi,h2, s. In principle, we could have followed the same 
procedure as in section ^ to obtain the dominant two loop corrections, i.e. 
replacing hi by h[ and ht by ht{mt). However, in order to obtain numerically 
correct results also in the regime Mgusy < 1 TeV, we did not expand V^^^ in 
powers of hi/Msusy, i-e. we used the full expression ( p.3|) for V^*^^). Then it 
becomes inconvenient to perform the field redefinition ( |3.6| ), which is implic- 
itly non- linear due to the hi dependence of t via nit. Therefore we proceed 
differently: For a given set of low energy parameters, which are implicitly 
obtained at the scale Q ~ Mgusy, we minimize directly 

V,jj = \/W + 1^(1) + (4.3) 
with as in (12.21), as in ( |2.3| ) and V^"^^ as in (|2.12| ). Points in the pa- 



rameter space leading to deeper unphysical minima of the effective potential 
with hi = or s = are removed. 

The overall scale is then fixed by relating the vevs hi to the physical Z 
mass through 

Ml = l{9! + 9l){Zj,hl + hl) (4.4) 

with Zh^ as in eq. ( p.l3| ). Next, we throw away all points in the parameter 
space where the top quark mass (including corrections ([2.14]) to ht) does 



not correspond to the measured mf"'*^ = 173.8 ± 5.2 GeV. We also ask for 
sfermions with masses mj > Mz/2 and gluinos with masses m~ > 200 GeV. 

Finally, the correct 3x3 Higgs mass matrix is related to the matrix of sec- 
ond derivatives of the Higgs potential at the minimum after dividing \d1^Veff 

by Zh^, and |f^/ii<9/i2 K// and ^dh^dgVeff by Z^^^. For each point in the param- 
eter space, we then obtain the two loop Higgs boson masses and couplings 
to gauge bosons. Then, we apply present constraints from negative Higgs 
search at LEP (cf. section ^for details). 

The results in section 5 are based on scannings over ~ 10^ points in the pa- 
rameter space. The essential effect of all constraints within the C(M-|-1)SSM 
is to further reduce the allowed range for the Yukawa coupling A to A < .3. 

4.2 General (M+1)SSM 

In the general (M+1)SSM, we only assume that we are in a local minimum 
of the effective potential ([4. 3]) and the running Yukawa couplings A, k, ht are 
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free of Landau singularities below the GUT scale. In order to scan the high 
dimensional parameter space ( [4.1|) of the general (M+1)SSM we proceed as 
follows: 

First, we use the three minimization equations of the full effective poten- 
tial ( [4 .31) with respect to hi, /12 and s in order to eliminate the parameters 
m^^, m'jj^ and m| in favour of the three Higgs vevs. Using the relation ( [4 .41) , 
we replace hi,h2 by tan/5 and Mz- Finally, eqs. ( p. 141) , (|2.15|) and (|2.16|) 
allow us to express in terms of m^"'^ and the other parameters. 

We are then left with six 'tree level' parameters A, k, Ax, A^, s, tan/5, and 
three parameters appearing only through the radiative corrections, which we 
choose as At, Mgusy and 5, as defined in eqs. (|2.6| ) and 



Requiring that the Yukawa couplings are free of Landau singularities be- 
low the GUT scale and using the renormalization group equations of the 
(M-l-l)SSM 1^, one obtains upper limits on X,K,ht at the susy scale. The 
latter turns into a lower bound on tan/5 depending mainly on mf"'*^ and Mg^sy 
As expected from eq. ( p.ll| ), we observe that upper limits on Higgs masses 
are obtained when A is maximal. From the renormalization group equations, 
one finds that the upper limit on A increases with decreasing k, thus we 
choose K ~ and A = Xmax ~ -7 (which still depends on ht, i.e. on tan/5). 

As already mentionned, one can see from eq. ( p.ll|) that the lightest Higgs 
mass is maximized for moderate values of tan/?. Hence, except in fig. 4 where 
tan/? varies, we fix tan/5 = 2.7 which, as we shall see, maximizes the Higgs 
masses for mf°'^ = 173.8 GeV. 

Unless stated otherwise, the upper limits on the Higgs masses presented 
in the next section are given in the maximal mixing scenario [At = ^/QMsusy)■ 
We have also found that Higgs masses are maximized for small values of 6 
and fixed 6 = (thus ttiq = rriT = Mgusy)- In order to obtain the results 
presented in the next section, we have used numerical routines to maximize 
the Higgs masses with respect to the remaining three parameters A\, A,^, s. 



5 Reduced couplings versus mass bounds 

Let us start with the mass mi of the lightest Higgs scalar, independently 
of its coupling to gauge bosons. The upper limit on mi in the general 
(M-l-l)SSM is plotted in fig. 1 (straight line) as a function of Msusy (for 
^poZe _ j^-j^g g QgV). This limit is well above the one of the MSSM because 
of the additional tree level contribution to mi proportional to A^M| (cf. 
eq. (|1.3|)). At Mg^sy = 1 TeV we have mi < 133.5 GeV (in agreement with 
the analytic approximation ( p.ll|) ); at Msusy = 3 TeV this upper limit in- 
creases only by ~ 3 GeV. This weak dependence on Mgusy is due to the 
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negative two loop contributions to mi. 

Within the C(M+1)SSM, the combined constraints on the parameter 
space require A to be small, A < .3 [0,0. Accordingly, the upper limit 
on mi is very close to the one of the MSSM. It is shown as crosses in fig. 1, 
and reaches 120 GeV at M^usy = 1 TeV. In the following, we shall assume 
M,usy = 1 TeV. 

A sideremark on the behavior for small Mgusy is in order: From eq. ( p. 71) , 
it is obvious that, in the assumed limit 5 — 0, the assumption of maximal 
stop mixing (At = V^susy) cannot be maintained for 

mt < Msusy < (5.1) 

because it would imply a negative stop mass squared. Therefore, in the 
general (M+1)SSM, we choose in this regime such that the lightest stop 
mass squared remains positive. On the other hand, within the C(M+1)SSM, 
where soft susy breaking terms are related, the limit Msusy small is not 
leasable since it would contradict the negative results on sparticle searches. 

As discussed in the introduction, the upper limit on mi is not necessarily 
physically relevant, since the coupling of the lightest Higgs to the Z boson 
can be very small. Actually, this phenomenon can also appear in the MSSM, 
if sin^(/5 — a) is small. However, the CP odd Higgs boson A is then necessarily 
light (m^ ~ nih < Mz at tree level), and the process Z — > hA can be used 
to cover this region of the parameter space in the MSSM. In the (M+1)SSM, 
a small gauge boson coupling of the lightest Higgs Si is usually related to a 
large gauge singlet component, in which case no (strongly coupled) light CP 
odd Higgs boson is available. Hence, Higgs searches in the (M+1)SSM have 
possibly to rely on the search for the second lightest Higgs scalar 5*2. 

Let us now define Ri as the square of the coupling ZZSi divided by the 
corresponding standard model Higgs coupling: 

Ri = {Sii sin/3 + Si2 cos/5)' (5.2) 

where Sn, Si2 are the Hi,H2 components of the CP even Higgs boson Si, 
respectively. Evidently, we have < Ri < 1 and unitarity implies 

3 

i=l 

Fortunately, as it was already mentionned, in the extreme case i?i — > 
the upper limit on m2 is the same as the above upper limit on mi. On the 
other hand, scenarios with, e.g., i?i ~ i?2 ~ 1/2 are possible. In the following 
we will discuss these situations in detail. 
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We are interested in upper limits on the two lightest CP even Higgs bosons 
5*1^2- These are obtained in the limit where the third Higgs, 5*3, is heavy and 
decouples, i.e. i?3 ~ (This is the equivalent of the so called decoupling limit 
in the MSSM: the upper bound on the lightest Higgs h is saturated when the 
second Higgs H is heavy and decouples). Hence, we have Ri + R2 — 1- 

In the regime Ri > 1/2 experiments will evidently first discover the light- 
est Higgs (with nil < 133.5 GeV for Msusy = 1 TeV). The 'worst case sce- 
nario' in this regime corresponds to mi ~ 133.5 GeV, Ri ~ 1/2; the presence 
of a Higgs boson with these properties has to be excluded in order to test 
this part of the parameter space of the general (M+1)SSM. 

The regime where Ri < 1/2 (and hence 1/2 < i?2 < 1) is more delicate: 
here the lightest Higgs may escape detection because of its small coupling, 
and it may be easier to detect the second lightest Higgs. In fig. 2 we show 
the upper limit on m2 as a function of R2 in the general (M-l-l)SSM as a 
thin straight line. For R2 1 (corresponding to Ri —>■ 0) we obtain the 
announced result: the upper limit on a Higgs boson with i? — > 1 is always 
given by the previous upper limit on mi, even if the corresponding Higgs 
boson is actually the second lightest one. The same applies, of course, to 
the C(M-|-1)SSM where the upper limit on m2 is also indicated as crosses in 
fig. 2. In the following we will discuss this 'delicate' regime, i?i < 1/2 and 
1/2 < i?2 < 1; in some detail: 

Fortunately, one finds that the upper limit on m2 is saturated only when 
the mass mi of the lightest Higgs boson tends to 0. Clearly, one has to 
take into account the constraints from Higgs boson searches which apply 
to reduced couplings R < 1/2 - i.e. lower limits on mi as a function of 
/2i ^ 1 — i?2 ^ in order to obtain realistic upper limits on m2 vs R2- 

Lower limits on mi as a function of Ri (in the regime Ri < 1/2) have 
been obtained at LEP We use the following analytic approximation for 
the constraints on Ri vs mi in this regime: 

The resulting upper limit on m2 is shown in fig. 2 as a thick straight line. This 
constraint is automatically included in the C(M+1)SSM results (crosses). 
Present and future Higgs searches at LEP will lead to more stringent con- 
straints in the regime .1 < Ri < 1/2 [|l^. We approximate the possible 



constraints from a run at 198 GeV cm. energy and 200 pb"^ by 

The resulting upper limit on m2 is shown in fig. 2 as a thick dashed line. 
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It would be desirable to have the upper limit on m2 in the general 
(M+1)SSM for arbitrary lower limits on mi as a function of Ri. To this 
end we have produced fig. 3. The different dotted curves show the upper 
limit on as a function of R2 for different lower limits on mi (as indicated 
on each curve) as a function of Ri (as indicated at the top of fig. 3). 

In practice, fig. 3 can be used to obtain upper limits on the mass m2, in the 
regime -Ri < 1/2, for arbitrary experimental lower limits on the mass mi: For 
each value of the coupling i?i, which corresponds to a vertical line in fig. 3, 
one has to find the point where this vertical line crosses the dotted curve 
associated to the corresponding experimental lower limit on mi. Joining 
these points by a curve leads to the upper limit on 1712 as a function of 
i?2- We have indicated again the present LEP limit ( |5.4| ), already shown 
in fig. 2, which excludes the shaded region (m2 > 172.5 GeV for R2 = .5, 
m2 > 150 GeV for R2 = .75, etc). We have also shown again the possible 
LEP2 constraints on m2 arising from ( p.5[ ) as a thick dashed line. 

Lower experimental limits on a Higgs boson with R > 1/2 restrict the 
allowed regime for m2 (for i?2 > 1/2) in fig. 3 from below. The present lower 
limits on m2 from LEP are not visible in fig. 3, since we have only shown the 
range m2 > 130 GeV. Possibly Higgs searches at the Run II of the Tevatron 
push the lower limits on m2 upwards into this range. This would be necessary 
if one aims at an exclusion of the 'delicate' regime of the (M+1)SSM: Then, 
lower limits on the mass m2 - for any value of R2 between 1/2 and 1 - of 
at least 133.5 GeV are required; the precise experimental lower limits on m2 
as a function of R2, which would be needed to this end, will depend on the 
achieved lower limits on mi as a function of -Ri in the regime -Ri < 1/2. 

In principle, from eq. ( |5.3| ), one could have R2 > Ri with R2 as small 
as 1/3. However, in the regime 1/3 < R2 < 1/2, the upper bound on m2 
as a function of -R2 for different fixed values of mi can only be saturated if 
Ri = R2- Then it is sufficient to look for a Higgs boson with a coupling 
1/3 < R < 1/2 and a mass m < 133.5 GeV to cover this region of the 
parameter space of the (M+1)SSM. 

Finally, we consider the dependence of the upper bounds on mi, 2 on tan/3 
and the top quark pole mass. In fig. 4 we plot the upper limit on mi, 2 
(for i?i,2 = 1) against tan/3 for mf"''^ = 173.8 GeV as a thick straight line. 
Remarkably, as announced before, this tan/? dependence is very different from 
the MSSM: the maximum is assumed for tan/3 ~ 2.7 (with mi, 2 — 133.5 GeV 
in agreement with figs. 2,3). The origin of this tan/? dependence is the tree 
level contribution ~ sin^/3 to (3.11). The height and the location of the 
maximum varies somewhat with mf°'^; the thick dashed and dotted curves 
correspond to mf°'^ = 173.8 ±5.2 GeV, respectively. The absolute maximum 
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is at tan/5 ~ 3 with mi_2 — 135 GeV. 

In the 'dehcate' regime, where one has to search for the second hghtest 
Higgs with i?2 between 1/2 and 1, one could worry whether the tan/5 depen- 
dence of the upper hmit on m2 is different. This is not the case: As a thin 
straight hne we show the upper hmit on m2 in the extreme case R2 = 1/2 and 
173.8 GeV (where the LEP constraint ( p.4|) is taken into account), 



pole 
t 



which assumes again its maximum for tan/5 ^ 2.7 (now with m2 — 172.5 GeV 
in agreement with figs. 2,3). As above, the thin dashed and dotted curves 
correspond to mf°''^ = 173.8 ± 5.2 GeV, respectively, and the absolute maxi- 
mum is at tan/3 ~ 3 with m2 ~ 175.5 GeV. Within the C(M-M)SSM, where 
A is small, the dependence of the upper limit on m2 on tan/3 ressembles more 
to the one of the MSSM as shown as crosses in fig. 4. 

To conclude, we have studied the CP even Higgs sector of the general 
(M-l-l)SSM and the C(M+1)SSM including the dominant two loop correc- 
tions to the effective potential. We have emphasized the need to search for 
Higgs bosons with reduced couplings, which are possible within this model. 
Our main results are presented in fig. 3, which allows to obtain the constraints 
on the Higgs sector of the model both from searches for Higgs bosons with 
weak coupling {R < 1/2), and strong coupling (i? > 1/2). The necessary 
(but not sufficient) condition for testing the complete parameter space of the 
(M+1)SSM is to rule out a CP even Higgs boson with a coupling 1/3 < i? < 1 
and a mass below 135 GeV. The sufficient condition (i.e. the precise upper 
bound on m2 vs R2) depends on the achieved lower bound on the mass of 
a 'weakly' coupled Higgs (with < i? < 1/2) and can be obtained from 
fig. 3. At the Tevatron this would probably require an integrated luminosity 
of up to 30 fb~^ |]T^. If this cannot be achieved, and no Higgs is discov- 
ered, we will have to wait for the results of the LHC in order to see whether 
supersymmetry beyond the MSSM is realized in nature. 
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Figure Captions 



Figure 1: Upper limits on the mass mi of the hghtest CP even Higgs bo- 
son versus Mgusy in the general (M+1)SSM (straight line) and the 
C(M+1)SSM (crosses). 

Figure 2: Upper limits on the mass m2 of the second lightest CP even Higgs 
(in the regime R2 > 1/2) against R2 in the general (M+1)SSM (thin 
straight line); the general (M-l-l)SSM with LEP constraints (|5.4|) (thick 
straight line); the general (M+1)SSM with expected LEP2 constraints 
( p3|) (thick dashed line); the C(M+1)SSM with LEP constraints ( ^) 
(crosses) . 

Figure 3: Upper limits on the mass m2 against R2, for different lower limits 
on the mass mi (as indicated on each line in GeV) of the lightest Higgs 
boson for 1/2 < R2 < 1. Ri = 1 — R2 is shown on the the top axis. 
The boundary of the shaded area corresponds to the thick line in fig. 
also the dashed line is the same as in fig. ^ 

Figure 4: Upper limits on mi 2 with R12 = 1 (thick lines), and upper limits 
on m2 with R2 = 1/2 (thin lines) versus tan (3 in the general (M+1)SSM 
for mf'' = 173.8 GeV (straight), 179 GeV (dashed) and 168.6 GeV 
(dotted); upper limit on m2 in the C(M+1)SSM (crosses) for mf"'"^ = 
173.8 ± 5.2 GeV. The LEP constraints ( |5.4|) are taken into account in 
each case. 
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